Search for: All records

The probability of finding a spherical “hole” of a given radius r contains crucial structural information about many-body systems. Such hole statistics, including the void conditional nearest-neighbor probability functions GV(r), have been well studied for hard-sphere fluids in d-dimensional Euclidean space Rd. However, little is known about these functions for hard-sphere crystals for values of r beyond the hard-sphere diameter, as large holes are extremely rare in crystal phases. To overcome these computational challenges, we introduce a biased-sampling scheme that accurately determines hole statistics for equilibrium hard spheres on ranges of r that far extend those that could be previously explored. We discover that GV(r) in crystal and hexatic states exhibits oscillations whose amplitudes increase rapidly with the packing fraction, which stands in contrast to GV(r) that monotonically increases with r for fluid states. The oscillations in GV(r) for 2D crystals are strongly correlated with the local orientational order metric in the vicinity of the holes, and variations in GV(r) for 3D states indicate a transition between tetrahedral and octahedral holes, demonstrating the power of GV(r) as a probe of local coordination geometry. To further study the statistics of interparticle spacing in hard-sphere systems, we compute the local packing fraction distribution f(ϕl) of Delaunay cells and find that, for d ≤ 3, the excess kurtosis of f(ϕl) switches sign at a certain transitional global packing fraction. Our accurate methods to access hole statistics in hard-sphere crystals at the challenging intermediate length scales reported here can be applied to understand the important problem of solvation and hydrophobicity in water at such length scales. 

We identify several distinct phases of thermalization that describe regimes of behavior in isolated, periodically driven (Floquet), mesoscopic quantum chaotic systems. In doing so, we also identify a Floquet thermal ensemble, the “ladder ensemble,” that is qualitatively distinct from the “featureless infinite-temperature” state that has long been assumed to be the appropriate maximum-entropy equilibrium ensemble for driven systems. The phases we find can be coarsely classified by (i) whether or not the system irreversibly exchanges energy of order ω with the drive, i.e., Floquet thermalizes, and (ii) the Floquet thermal ensemble describing the final equilibrium in systems that do Floquet thermalize. These phases are representative of regimes of behavior in mesoscopic systems, but they are sharply defined in a particular large-system limit where the drive frequency ω scales up with system size N as the N → ∞ limit is taken: we examine frequency scalings ranging from a weakly N-dependent ω(N)∼logN, to stronger scalings ranging from ω(N)∼√N to ω(N)∼N. We show that the transition where Floquet thermalization breaks down happens at an extensive drive frequency and, beyond that, systems that do not Floquet thermalize are distinguished based on the presence or absence of rare resonances across Floquet zones. We produce a thermalization phase diagram that is relevant for numerical studies of Floquet systems and experimental studies on small-scale quantum simulators, both of which lack a clean separation of scales between N and ω. A striking prediction of our work is that, under the assumption of perfect isolation, certain realistic quench protocols from simple pure initial states can show Floquet thermalization to a type of Schrodinger-cat state that is a global superposition of states at distinct temperatures. Our work extends and organizes the theory of Floquet thermalization, heating, and equilibrium into the setting of mesoscopic quantum systems. 

The out-of-time-ordered correlator (OTOC) has emerged as an interesting object in both classical and quantum systems for probing the spatial spread and temporal growth of initially local perturbations in spatially extended chaotic systems. Here, we study the (classical) OTOC and its “light cone” in the nonlinear Kuramoto-Sivashinsky (KS) equation, using extensive numerical simulations. We also show that the linearized KS equation exhibits a qualitatively similar OTOC and light cone, which can be understood via a saddle-point analysis of the linearly unstable modes. Given the deep connection between the KS (deterministic) and the Kardar-Parisi-Zhang (KPZ, which is stochastic) equations, we also explore the OTOC in the KPZ equation. While our numerical results in the KS case are expected to hold in the continuum limit, for the KPZ case it is valid in a discretized version of the KPZ equation. More broadly, our work unravels the intrinsic interplay between noise/instability, nonlinearity, and dissipation in partial differential equations (deterministic or stochastic) through the lens of OTOC. 

We investigate ergodicity, chaos, and thermalization for a one-dimensional classical gas of hard rods confined to an external quadratic or quartic trap, which breaks microscopic integrability. To quantify the strength of chaos in this system, we compute its maximal Lyapunov exponent numerically. The approach to thermal equilibrium is studied by considering the time evolution of particle position and velocity distributions and comparing the late-time profiles with the Gibbs state. Remarkably, we find that quadratically trapped hard rods are highly nonergodic and do not resemble a Gibbs state even at extremely long times, despite compelling evidence of chaos for four or more rods. On the other hand, our numerical results reveal that hard rods in a quartic trap exhibit both chaos and thermalization, and equilibrate to a Gibbs state as expected for a nonintegrable many-body system. 

We study one-dimensional hybrid quantum circuits perturbed by quenched quasiperiodic (QP) modulations across the measurement-induced phase transition (MIPT). Considering non-Pisot QP structures, characterized by unbounded fluctuations, allows us to tune the wandering exponent β to exceed the Luck bound ν ≥ 1/(1−β) for the stability of the MIPT, where ν = 1.28(2). Via robust numerical simulations of random Clifford circuits interleaved with local projective measurements, we find that sufficiently large QP structural fluctuations destabilize the MIPT and induce a flow to a broad family of critical dynamical phase transitions of the infinite QP type that is governed by the wandering exponent β. We numerically determine the associated critical properties, including the correlation length exponent consistent with saturating the Luck bound, and a universal activated dynamical scaling with activation exponent ψ ≅ β, finding excellent agreement with the conclusions of real-space renormalization group calculations. 

Random quantum circuits continue to inspire a wide range of applications in quantum information science and many-body quantum physics, while remaining analytically tractable through probabilistic methods. Motivated by an interest in deterministic circuits with similar applications, we construct classes of nonrandom unitary Clifford circuits by imposing translation invariance in both time and space. Further imposing dual unitarity, our circuits effectively become crystalline spacetime lattices whose vertices are swap or iswap two-qubit gates and whose edges may contain one-qubit gates. One can then require invariance under (subgroups of) the crystal’s point group. Working on the square and kagome lattices, we use the formalism of Clifford quantum cellular automata to describe operator spreading, entanglement generation, and recurrence times of these circuits. A full classification on the square lattice reveals, of particular interest, a “nonfractal good scrambling class” with dense operator spreading that generates codes with linear contiguous code distance and high performance under erasure errors at the end of the circuit. We also break unitarity by adding spacetime translation-invariant measurements and find a class of such circuits with fractal dynamics. 

Many-body localized (MBL) systems fail to reach thermal equilibrium under their own dynamics, even though they are interacting, nonintegrable, and in an extensively excited state. One instability toward thermalization of MBL systems is the so-called “avalanche,” where a locally thermalizing rare region is able to spread thermalization through the full system. The spreading of the avalanche may be modeled and numerically studied in finite one-dimensional MBL systems by weakly coupling an infinite-temperature bath to one end of the system. We find that the avalanche spreads primarily via strong many-body resonances between rare near-resonant eigenstates of the closed system. Thus we find and explore a detailed connection between many-body resonances and avalanches in MBL systems.